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Abstract 

We consider a lattice version of the bosonized Gross-Neveu model. 
It is explicitely chiral symmetric and its numerical simulation does 
not involve any anticommuting field. We study its non trivial 1/N 
expansion up to the next-todeading term comparing the results with 
explicit numerical simulations. 
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1. Introduction 

It is a well known property of two dimensional space-time that a free 
massless fermion is equivalent to a free massless boson fl[]. At the operator 
level, we can find a mapping, the Mandelstam representation [Q], which solves 
the fermionic field equations in terms of bosonic operators. The vector and 
axial currents, together with the energy momentum tensor satisfy the same 
Kac-Moody and Virasoro algebras in the two formulations || allowing for 
a full correspondence. The same renormalized Green functions are obtained 
for the usual bilinear local fermionic operators and suitable local bosonic op- 
erators. The mapping can be extended to interacting models. Bosonization 
is somewhat peculiar of 1 + 1 dimensions where the spin assignment is arbi- 
trary but, on the other hand, these models can help in studying physically 
interesting properties of realistic theories like confinement, chiral symmetry 
breaking and the baryonic spectrum. From the point of view of the non 
perturbative analysis, we show in this letter how the bosonization technique 
is useful in solving the problems encountered in the simulation of fermionic 
models. The expensive inversion of the fermionic propagator is a serious ob- 
stacle for numerical methods. Moreover, for a large class of lattice models, 
the interplay between fermion doubling and chiral symmetry [ffl complicates 
the analysis of the continuum limit. The bosonization approach is, in our 
opinion, powerful in both respects. It avoids Grassmann fields ab initio and 
may potentially overcome the Nielsen- Ninomiya theorem as shown in |J . 
There are also additional advantages concerning the so-called Symanzik im- 
provement, but we postpone their discussion. In this letter, we propose a 
lattice version of the two dimensional Gross-Neveu model. It is a favoured 
theoretical laboratory for the analysis and simulation of fermionic models be- 
cause it has nice analytical properties: it dynamically breaks a discrete chiral 
symmetry and is 1/N expandable allowing for an important check of the sim- 
ulation m [7|, |J. Its spectrum of bound states is described by a semiclassical 
formula from which mass ratios can be extracted. We shall not consider here 
the problem of non-Abelian bosonization || and the number of flavors N 
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will be only an expansion parameter. As a consequence, the largest part of 
the model's symmetry will be non linearly realized on the bosonic fields. 



2. The continuum model 

The 1 + 1 free massive fermion action may be bosonized in terms of the 
Mandelstam representation which expresses left and right handed Weyl 
fermions using a scalar massless field (p(x): 

^ = v^ :exp { _ ^(£c de7r(e)+ ^ (x) )} : ' (1) 

= v^ :exp {"^(roo de7r(e) "^ (:c) )} : • (2) 

where n(x) = d (p(x) and c = |e 7 is related to the short distance expansion 
of the two dimensional free bosonic propagator. The normal ordering flU 
performed with respect to the mass \x []. The bosonization dictionary is 



is 



#^ -»• n(9a<f) 2 — :cos/fy: (3 2 = Att. (4) 

2 7T 

Actually, the above construction can be extended to the massive Thirring 
model and Eqs. [| correspond to the limit of vanishing interaction. In the 
general case, (3 is a function of the four-fermion coupling. Our starting ob- 
servation is that the purely fermionic Gross- Neveu model may be formulated 
recasting the four fermion interaction in the form of a space-time dependent 
mass term. This suggests a possible bosonization procedure which is partic- 
ularly useful for an expansion in the number of flavours. In the continuum, 
the fermionic Gross-Neveu model is described by the tree level action 

S = J d 2 x (V Q) <^ (a) - ^g 2 (^ a V (a) ) 2 ) a = l---JV, (5) 



We recall that in a continuum model with the cut-off scale A, we have 

: cos Pip := — cos (Hip (3) 
A* 

where fi is an arbitrary renormalization scale. 



2 



where ip^ a ' is a multiplet of two-dimensional spinors. Rewriting the quartic 
interaction with a Lagrange multiplier (see []TTJ for a complete discussion of 
the model renormalization in this formulation) we obtain 

S = Jd 2 X (V^W^ + <T$ (a ty (a) + • ( 6 ) 

The above dictionary for bosonization can be used. At the level of formal 
manipulations, the a field is just a combinatorial device and our scheme is 



equivalent to Witten's treatment of the bosonized Gross-Neveu model |L2 
|13| , |I4}| . We now propose a lattice model and show several indications that 
it can be used to study the bosonic sector of the Gross-Neveu model, but it 
must be kept in mind that the precise nature and existence of its continuum 
limit is an open problem. 

3. The lattice model 

We want to follow the analysis in ||. Therefore, we write on the lattice 

5 = E^n + Sa$) + Sint(a,$) , (7) 

where $ is a set of fields which we now specify. For the fermionic model 
of 0, $ is just a multiplet of iV Dirac fermions ip( a \ S^ n (-ip) is some lattice 
kinetic action which avoids doubling and 

S- mt (a^) = Yl CC^n- (8) 

n,ce 

For the bosonized model, $ is a set of N massless bosonic fields Ski n (v) 
is 

SkinM = I A+< A+< A+< = < +A - < (9) 

and 

S in t{cr,<t>) = -u^a n cos/3y3° . (10) 
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The constant 10 comes from a tadpole computation. In a free scalar theory 
with mass fi one can write 



<cos/^ n ) = e-^ A «." , (11) 

and therefore 

1 (3 2 f d 2 p 1 

H : cos (3<p : = -R{fi) cos /3(p R(ji) = /iexp — / — - . (12) 

CI — •> \Z/J^ 1 — p ^Li 

We take 



u = = ^— ~ 1.6035286 . (13) 

7T 7T 

The symmetry 

7T 

a^-a p^ip + - (14) 

is unbroken at the level of Feynman rules and corresponds to the following 
symmetry of the fermionic realization 

%j) -> 75^ ip -> -^75 cr -> -a. (15) 

This symmetry is dynamically broken: the cr field acquires a non zero expec- 
tation value which explains infrared stability. As in the fermionic model, 
volume effects may be strongly reduced utilizing the Symanzik improve- 
ment fT5 l- Here comes another advantage of the bosonized model: the 



Symanzik improvement adds no more difficulties to the numerical simulation. 
On the other hand, in the fermionic model, the additional interactions reduce 
the extent to which the fermionic matrix is sparse. As a consequence, pre- 
conditioning becomes extremely more difficult. As an aside, we remark that 
an alternative approach to the introduction of the auxiliary field would be to 
bosonize the purely fermionic model. The operator product expansion of the 
four-fermion interaction alters the kinetical term and fixes a trajectory in the 
(A, (3) plane along which the critical point is reached |16| . In |1^, [13], [14]], this 



kind of bosonization has been applied to the Gross-Neveu model in order to 
obtain important analytical results in the continuum at particular values of 
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N. Here, we have chosen to introduce the field a to mimic the 1/N expansion 
of the purely fermionic model. 



4. The N = 00 limit 

At leading order, the effective action for the a field is 

r((7) 4"^ 0g(e " SW> (16) 

where L is the lattice side, a is kept space-time independent and the vacuum 
expectation value is computed with respect to the free bosonic action. In the 
fermionic model, the interaction action is simply a mass term and therefore 
the second term in T can be explicitely evaluated 

1 1 

— log(e- 5 -) = — log det(A^ + a) = f log(A^ + a) . (17) 
If we have instead 

Sint = COS /Vn (18) 

n 

then a straightforward computation is not possible any more. However, we 
can consider 

F eff (a) = -r (a) = F d (<r) + F qu (a) (19) 
Fci(a) = ~ (20) 

Fqu(o-) = -^(Ecosfc) (21) 

where the vacuum expectation value is taken now with respect to the action 
of the N — 1 model. A Monte Carlo calculation can give very precise results 
without being too expensive. We expect to find 

Urn F qu (a) = ±uj F qu (a) = -F qu (-<r) . (22) 

<x — >±oo 

On a finite lattice F qu (o") has a finite slope as o — > which signals the 
existence of a finite volume critical coupling for chiral symmetry breaking. 
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Asymptotic scaling depends on the non trivial infinite volume, small a be- 
haviour of F qu (cr). We want to show that the solution of F e s(a) = satisfies 
the scaling law 

a = A B exp{-j) (23) 
This is true if we prove that the expansion 

F qu (cr) = const a - -a log a + 0(a 2 \og k a) (24) 

7T 

does hold. However, a naive expansion of (cos (3(p) in powers of a fails when 
L — > oo. The first terms are indeed 

(cos 0<p) = e~^ 2j ° +ujaJ2e' 02Jo (cosh(3 2 J n -l) + --- 

n 

= co(L) + oua Ci(L) H (25) 

where 

J„ = ^E^"4 ( 26 ) 

computed introducing antiperiodic boundary conditions. At a given a we 
cannot send L — > oo because Ci(L) — > oo. The question is: how can we 
relate the IR leading logarithms of the <r = 0, L — > oo limit to those of 
the a — * 0, L = oo one ? Let us consider a function F(cr, L) such that the 
F(a, oo) exists and is non analytic at a = 0. At fixed L, we can expand in 
powers of a 

F(a,L) = J2c a (L)a a (27) 

with vanishing radius of convergence as L — > oo. In this model, the volume 
effects can be kept constant uniformly in a by fixing the adimensional variable 
S = aL. Therefore the leading IR logs are buried in the expansion 

G(a) = J2c a (-) * a • (28) 
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This simple approach works if the arbitrary constant k is irrelevant for the 
determination of the logarithms. This is the case for F qu (a) []. A numerical 
investigation gives as L — ► oo 

/r . const 
Co{L) ~ — , (32) 

ci(L) ~ z/logL . (33) 

With v ~ 0.12379219 from which it follows that u 2 v is equivalent to 1/tt 
up to our precision. Besides this strong numerical indication of asymptotic 
scaling in the N = oo limit, we give a sketchy analytical confirmation that 
the slope is indeed the correct one. Let us start from 

X = 2e-^£smh 2 ^J„. (34) 

n * 

If we introduce a regularization mass /i 2 and go to the continuum in the 
second factor 

X 2 exp — /? 2 J T^^I" / rf 2 xsinh 2 ^(0(^(0)) . (35) 
The leading singularity in the // — > limit is 

(0(x)0(O)) ~ — — log Vx^cfx (36) 

27T 



2 Another example is the case of cq(L). A numerical evaluation gives with great accu- 
racy 

, . const 

C(L) ~ — . (29) 

We can assume 

J2 1~ const + i- log L. (30) 

p 

Then, if one substitutes L ~ the correct logarithm in the exact expansion of 



d 2 p 1 



(2tt) 2 j5 2 + M 2 



(31) 



as a function of p is indeed obtained. 
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whence, formally 



X = 



J £x h ■ (37) 



2tt 2 uj 2 

We integrate this ill-defined quantity between the circles p = a and p = L to 
pick up the volume dependence and get 

X = — 2 \og-. (38) 



5. 1/N corrections 

The next to leading corrections in the 1/N expansion take the form 

r(<0 = ^-^;log(e- s «)+ (39) 



2N J (2tt) 

n(p, a) is the one particle irreducible a self energy computed with the N — 1 
action in which the dynamical field a is replaced by its N = oo value. In the 
fermionic model, II (p, a) is just a one loop lattice Feynman diagram. In the 
bosonized model, the numerical determination of U(p, a) is not straightfor- 
ward. We must utilize a Monte Carlo simulation to obtain the a field two 
point function. 

- cj 2 ((cos f3ip(x) cos /3(p(y)) - (cos fi^) 2 ) (40) 

giving just U(x — y,a). By Fourier transforming, we obtain U(p, a) and can 
compute the desired 0(1/ N) correction. We then write 

r(*) = ^ + rx((7) + ^r 2 (a) + o(^) (4i) 

and expand in powers of 1/iV the solution of r'(cr) = 0: 

^ = ^o + ^i + 0(^). (42) 
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The desired correction is therefore 

° l = -i/a-f;,M • < 43 > 

It is odd in a consistently with the unbroken discrete chiral symmetry. 

6. The simulation 

The simulation of the model has been performed using the Hybrid Monte 
Carlo algorithm |17]]. As a first step, we have determined F qu (a) on a 40 2 



lattice varying a over a wide range of values using the N — 1 model. The 
resulting plot is shown in Fig. 1. At L = 40, the critical coupling for chiral 
symmetry breaking is about A c ~ 0.38. Using the couplings A which can 
be extracted from that figure, we have studied the behaviour of F qu (a) and 
r 2 (cj) around o = 0.2, 0.25, 0.3 in order to determine with sufficient numerical 
accuracy the quantities needed to compute o\. Then, we have simulated the 
full model with iV = 20 flavours. In Tab. [I], we show at three values of A, 
the leading 1/N prediction for (a), its 1/N correction, the 1/N corrected 
value considering N = 20 flavors and the explicit numerical data including 
the measure of (cr 2 ) c = (cr 2 ) — (cr) 2 . We see that the leading correction 
is necessary. The theoretical error on the correction is rather large; the 
discrepancy with the numerical data is therefore about (8 ± 5)%. All these 
simulations has been performed using the 40 2 lattice and a non-improved 
action. We also simulated at the lower A = 0.49 utilizing a 60 2 lattice. We 
did not compute the theoretical 1/N corrections, but used this point just to 
look at the behaviour of 

i?,(A) = (a)exp(-^-^) (44) 



N-1X 

which is also reported in Tab. [I]. Having at our disposition two different 
lattice models, we can relate the approximate finite renormalization shift in 
the coupling to the A parameters ratio. In 0, we find the value of A at 
which the Symanzik action gives (cr) = 0.1488, 0.2, 0.25 and 0.3 at N = oo, 
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L = oo. Using 

1 1 1 , A B 

)T-)T = ^ log / (45) 

A B Af ft Af 

and the result on the N = oo, L = AO lattice for the bosonized model, we 
obtain A B = 90.6, 85.0, 80.7, 75.2. The different results at different A are 
due to finite size effects at L = 40 (60 for the first point) with the non- 
improved action and to non universal terms of the lattice f3 function. The 
values \R — A B \/A B , around 10%, may be explained with the 1/N corrections 
affecting A B . 

7. Conclusions 

We have proposed a possible lattice action for a bosonization of the Gross- 
Neveu model. There are many advantages with respect to the fermionic 
formulation from the point of view of the Monte Carlo simulation. There 
are no anticommuting fields and fermion doubling can be naturally avoided. 
The computer codes are extremely fast and there is no need for any fine 
tuning in order to restore the Z2 chiral symmetry. In this letter we were 
not interested in studying scaling properties or the mass spectrum of the 
model. Instead, we have shown the great simplicity and efficiency of its 
numerical simulation including its 1/N expansion. A detailed description 
of the simulation including scaling results, a study of the lightest bosonic 
bound states and a thorough analysis of the Symanzik improved action will 
be given elsewhere. To be fair, we must remark once again that we limit our 
analysis to the bosonic sector. In our opinion, work on lattice actions for the 
bosonized QCD 2 could be extremely interesting. 
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Table 1: Summary of the results 



A 


O"0 


o\ 


(J + Oi/20 


<*> 




R 


0.49000 


0.1488 






0.090(2) 


0.0606(3) 


77(2) 


0.51905 


0.2000 


-1.7(2) 


0.12(1) 


0.127(1) 


0.0681(6) 


74.6(6) 


0.54376 


0.2500 


-2.1(2) 


0.15(1) 


0.166(1) 


0.0744(6) 


72.7(4) 


0.56876 


0.3000 


-2.4(2) 


0.18(1) 


0.209(1) 


0.0807(7) 


70.0(3) 



Figure 1: F qu (a) on the 40 2 lattice 
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